Path Integral Quantization and 
Riemannian-Symplectic Manifolds 

Sergei V. SHABANOV a '[] and John R. KLAUDER b 

a Institute for Theoretical Physics, FU-Berlin, Arnimallee 14, Berlin, 

D- 14195, Germany 

b Departments of Physics and Mathematics, University of Florida, Gainesville, 
Q\ • FL-32611, USA 

G\ 

Abstract 

We develop a mathematically well-defined path integral formalism for general 
symplectic manifolds. We argue that in order to make a path integral quantization 
covariant under general coordinate transformations on the phase space and involve 
a genuine functional measure that is both finite and countably additive, the phase 
space manifold should be equipped with a Riemannian structure (metric). A suitable 
method to calculate the metric is also proposed. 
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1. Most quantization formulations lead to correct results only when the system in ques- 
tion is expressed in one of a family of special (Cartesian) coordinate systems jj]]. However, 
for Hamiltonian systems on a generic symplectic manifold, the canonical variables do not 
admit such special coordinates. Moreover classical Hamiltonian dynamics exhibits covari- 
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ance under general coordinate transformations on the symplectic manifold, and, in this 
sense, is coordinate-free. This covariance is lost upon a formal phase-space path integral 
quantization. The reason is that the Hamiltonian action involves terms linear in the time 
derivatives and, therefore, its exponential (even in Euclidean time) does not induce a 
proper (normalizable, cr-additive) measure on the path space which is covariant under 
general coordinate transformations. Since the discovery of the Hamiltonian path integral 
by Feynman, the indisputable fact of its coordinate dependence has generally complicated 
any straightforward use of Hamiltonian path integrals in applications. Indeed, to justify 
any formal manipulation regarding the path integral measure, the operator formalism has 
generally been invoked. The aim of this letter is to resolve this long-standing problem and 
to establish a rigorous, coordinate-free path integral formalism for general Hamiltonian 
systems. 

In P| it was shown that the formal phase-space path integral measure can be made 
covariant under general coordinate transformations by means of a special regularization 
that involves an auxiliary Brownian motion on a flat phase space. The desired quantum 
mechanics is restored by taking the diffusion constant to infinity. The modified Hamilto- 
nian path integral determines a stochastic process on the phase space in which the original 
first-order Lagrangian plays the role of the external potential. The stochastic process on a 
flat phase space remains flat under coordinate transformations. Thus one may change the 
variables in the regularized path integral according to the rules of the stochastic integral 
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calculus, and at the end of calculations one has to take the infinite diffusion constant 
limit. 

To be specific, and also to introduce notation, consider quantum mechanics in a flat 
phase space. Let g J be a set of Cartesian coordinate canonical operators satisfying the 
Heisenberg commutation relations [q\q k ] = i uj ■ ?fe . Here uj ^ k — — uj fcj is the canonical 

. o , 

symplectic structure. We introduce the canonical coherent states as \q) = |0), 
where uj j n uj nk = 6 k , and |0) is the ground state of a harmonic oscillator with unit angular 
frequency. Any state |?/>) is given as a function on phase space in this representation by 
(q\ip) = ip(q)- A general operator A can be represented in the form A = J dqa(q)\q)(q\, 
where a(q) is the lower symbol of the operator and dq is a properly normalized form of 
the Liouville measure. The function A(q, q') = (q\A\q') is the kernel of the operator. 

The main object of the path integral formalism is the integral kernel of the evolution 
operator 

<?(*)=<? 



K t (q,q') = (q\e- UH \q') 



Vqe 



J*drfiq^ jk q k -hj 



q(0)=q' 



Here H is the Hamiltonian, and h(q) its symbol. The measure formally implies a sum 
over all phase-space paths pinned at the initial and final points, and a Wiener measure 
regularization implies the following replacement 



Vq^Vfi u (q) =Vqe-^ti dT * 2 = N u (t) dfx v w (q) 



(2) 



The factor N v (t) equals 2ne ut/2 for every degree of freedom, d^ v w [<q) stands for the Wiener 
measure, and v denotes the diffusion constant. We denote by K"(q, q') the integral kernel 
of the evolution operator for a finite v. The Wiener measure determines a stochastic 
process on the flat phase space. The integral of the symplectic one-form / q uj dq is a 
stochastic integral that is interpreted in the Stratonovich sense. Under general coordinate 
transformations q = q(q), the Wiener measure describes the same stochastic process on 
flat space in the curvilinear coordinates dq 2 = da(q) 2 , so that the value of the integral is 
not changed apart from a possible phase term. After the calculation of the integral, the 
evolution operator kernel is obtained by taking the limit v — ► oo. The existence of this 
limit, and also the covariance under general phase-space coordinate transformations, can 
be most easily proved through the operator formalism for the regularized kernel K"(q, q'). 

Note that the integral ([!]) with the Wiener measure inserted can be regarded as an 
ordinary Lagrangian path integral with a complex action, where the configuration space 
is the original phase space and the Hamiltonian h(q) serves as a potential. Making use of 
this observation it is not hard to derive the corresponding Schrodinger-like equation 



UJ 



■<kQ 



subject to the initial condition K^ =0 (q,q') = 8(q 



ih(q) 



q'),0 



K?(q,q') 



(3) 



< v < oo. One can easily 



show that K" — > K t as v — ► oo for all t > 0. The covariance under general coordinate 



2 



transformations follows from the covariance of the "kinetic" energy of the Schrodinger 
operator in (||): The Laplace operator is replaced by the Laplace-Beltrami operator in the 
new curvilinear coordinates q = q(q), so the solution is not changed, but written in the new 
coordinates. This is similar to the covariance of the ordinary Schrodinger equation and 
the corresponding Lagrangian path integral relative to general coordinate transformations 
on the configuration space: The kinetic energy operator (the Laplace operator) in the 
ordinary Schrodinger equation provides a term quadratic in time derivatives in the path 
integral measure which is sufficient for the general coordinate covariance ||. We remark 
that the regularization procedure based on the modified Schrodinger equation (|3]) applies 
to far more general Hamiltonians than those quadratic in canonical momenta and leading 
to the conventional Lagrangian path integral. 

The regularization procedure developed above deserves a further mathematical elab- 
oration. To this end we point out that the operator K\ acts in a Hilbert space, denoted 
here and belows as 7i u , that is larger than the original Hilbert space TCoo- Indeed, con- 
sider canonical pairs [Q k ,P n ] = i5 k , [Pk,P n ] = [Q k ,Q n ] = and sharp-eigenvalue states 
\q) (observe the difference in notation!) being eigenstates of the complete set of commu- 
tative operators: Q k \q) = q k \q)- These states form a basis in 7i v . So a resolution of unity 
reads \ v = J dq\q)(q\ where (q\q') = S(q — q'). In constrast to |g), the coherent states \q) 
form a basis in the smaller space 7ioo, though they carry the same label q. Note that 
the coherent states in Ti v would have twice as many labels, say p and q, associated with 
the canonical operators P and Q, respectively. In 7i v the original canonical operators 
can be represented as q k = — u kn P n + \Q k , that is, we have a reducible representation 
of the original Heisenberg algebra. It is also important to observe that the extension of 
any opearator in Ti^ to the larger space H u is by no means unique. This issue is not, 
however, our main concern. 

To achieve a continuous-time regularization of the path integral measure, and thereby 
to make the latter mathematically well defined, it is sufficient to consider a specific ex- 
tension of the unitary operator K t = exp(-itH) which is induced by the extension of the 
Hamiltonian (cf. (||)) 

H^H, = -^(q-Q) 2 + h(Q) , (4) 

where the function h is the (lower) symbol of the original Hamiltonian H in the canonical 
coherent state representation. Observe that h(Q) is uniquely defined, no operator ordering 
ambiguities occur since the operators Q k are commutative. In the operator formalism 
developed, the Wiener measure regularization is based on the important relation 

Jim(g|e-^|g') = Jim K?(q,q') = (q\e- itA \q') , t > . (5) 

Thus, the ill-defined formal path integral representation of the r.-h.s. is replaced by the 
convergent sequence of well-defined stochastic integrals in the l.-h.s. which, as has been 
argued above, are explicitely covariant under general coordinate transformations q k — > 
q k (q). Note also that Hi ^ H u , and the immaginary part of the extended Hamiltonian 
H u gives rise to the Wiener measure in the path integral representation. 
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Now let q l be local coordinates on a symplectic manifold, with 

{q^q k } q = ^\q) . (6) 

The symplectic two- form u> = u^dq 1 A dq 3 ', u>ikUJ kj = 5j, is closed duj = (the Jacobi 
identity for the Poisson bracket ([])). A formal path integral quantization ([I]) (where 
the symplectic volume element and the symplectic one-form are modified according to 
the symplectic structure (|5])) runs into notorious difficulties if one tries to relate it to 
the operator formalism. In particular, by means of the formal change of the integration 
variables to Darboux variables, the overlap function (q\q') (the integral kernel of the 
unit operator, h = 0) seems always possible to be expressed as the integral (p]) with the 
canonical symplectic structure. This is certainly not true. Coherent states on general 
manifolds do not coincide with the canonical coherent states. Therefore, the path integral 
measure must be regularized relative to general coordinate transformations on the phase 
space. A possible way to achieve this is to replace the formal, ill-defined Liouville measure 
by the Wiener measure with an appropriate metric. However, there is a great deal of 
freedom in choosing such a metric 0. The overlap function (q\q') and, hence, the entire 
Hilbert space will depend on this choice. 

The main goal of our work is to give a procedure of how the metric can be constructed 
and to show that the corresponding path integral is coordinate-free and determines a 
consistent quantum mechanics. 

2. To avoid the difficult problem of quantizing general symplectic structure (§), we 
shall use the conversion formalism |4], ||. The idea is to embed the symplectic manifold 
into a larger Euclidean (flat) phase space and impose special second-class constraints such 
that, when reduced on the constraint surface, the canonical symplectic structure of the 
embedding space turns into the original symplectic structure. The second-class constraint 
system is then replaced by an effective Abelian gauge theory, and the latter is quantized 
by the standard methods. The algebra of gauge invariant operators is then taken as 
a quantum algebra of observables associated with the symplectic structure (|j). The 
representation (Hilbert) space of this algebra is a subspace of gauge invariant functions 
in the effective quantum gauge theory. 

Let us turn to the details. Consider an extended Euclidean phase space spanned by 
the Cartesian coordinates q l and momenta pi, and equipped with a canonical symplectic 
structure which we choose as {q\q 3 } = {Pi,Pj} = and {q l ,Pj} = 6). The second-class 
constraints are 

(pj =p j + a j (q) = , (7) 

where ai(q) are coefficients of the symplectic one-form a = atdq 1 and da — u being 
the symplectic two-form. It is not hard to verify that Ay = {Lpi,(pj} = u)ij(q). The 
symplectic structure on the constraint surface is induced by the Dirac bracket {■, -}d = 
{-,-} — {-, ipi}A lJ {(pj, ■}, where A* fe A^- = <5j. Eliminating the variables Pi by means 
of the constraints and assuming the variables q l to be the coordinates on the physical 
phase space, we find that the induced symplectic structure coincides with the original 
one: {q J ,q k }D\ip=o = u> jk (q). Thus, a Hamiltonian dynamics on a generic symplectic 
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manifold can be described as a second-class constraint system with the action S v = 
I (PsO* ~ M<?) — -^Vj) dt on an extended phase space. 

To quantize the theory, the canonical variables pj and q 1 are promoted into self- 
adjoint operators satisfying the canonical commutation relations ([,]=&{,})■ One next 
introduces a new set of canonical variables Qi associated with the second-class constraints 
ifj also satisfying the Heisenberg commutation relations , 6 ] = % u J and commuting 
with <y? and p k . A general operator A can be represented in the form 

A = J2 A n-u& <?) s y m { §n ■ ■ ■ #0 = E 4%»& 9)%"'" > (8) 

n n 

where the symbol sym means a symmetrized product of the operators in the brackets. 
The physical subspace in the extended Hilbert space is selected by the gauge invariance 
condition 

^■|'</W = 0, (9) 

where the generators of gauge transformations &j and the Hamiltonian in the extended 
Hilbert space are constructed to satisfy the conversion equations j|, |J 

[a J ,a k ] = 0, [H g ,a J } = 0, (10) 

subject to the boundary conditions af^ = ipj, H^> = H(q). In the classical limit, the 
dynamics of the physical degrees of freedom of the so constructed gauge theory coincides 
with that on the original symplectic manifold ||, that is, the correspondence principle is 
fulfilled. 

To solve the conversion equations, one sets 0, 0] Gj = 0j + Aj(q, 6), A^\q) = 0. The 
operators g- 7 commute, so we omit the hat while analyzing the equation for Aj = Aj(q, 6) 
that results from (|10|) 

djA k - d k Aj - i[Aj, A k \ = uj jk (q) . (11) 

Here dj = d/dqK The solution is not unique. If A k is a solution then for any A = A(q, 8), 
the operator A£ = e lA A k e~ 1 ^ — ie l!v d k e~ l!v is also a solution, provided this transformation 
does not violate the zero boundary condition A k ^ = 0. In fact, this condition imposes 
some restrictions on A discussed later. 

An important point to observe is that the physical amplitudes do not depend on A. 
Indeed, any physical amplitude reads {ipi\0\ip 2 ) where ^1,2) are gauge invariant (physical) 
states determined by @ and O is any operator satisfying the condition [O, <7j] ~ <7j. If A k 
is replaced by A k , then the new Abelian constraints are a k = e lK a k e~ tk . Thus, the new 
sets of physical operators and states are obtained from the former ones by the unitary 
transformation e* A , that is, the physical amplitudes are not changed. It is therefore 
sufficient to solve the conversion equations modulo the A-transformations. 



3. The right-hand side of Eq. ((TTJ) can be regarded as the field strength in the Yang- 
Mills theory with an infinite dimensional gauge group generated by the operators e( n ). 
Introducing the Yang-Mills potential C k = a k + A k , we observe that the corresponding 
field strength vanishes if A k satisfies (0). This means that the vector potential C k must 
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be a pure gauge, C k = —ie in dke lQ , where Q = Q(q,9). Yet, this solution is determined 
modulo the gauge (A-)transformations. The operator Cl satisfies only the condition that 
C k = ctfc and this condition should not be violated under the gauge transformations 
of Ck- For this reason, the pure gauge potential Ck can not be gauged to a vanishing 
potential (which does not satisfy the condition C k = ct^). 

The operators e(o) and are generators of the Heisenberg algebra which is a finite 
dimensional subalgebra of the algebra generated by all the e( n ). So we can always represent 
the gauge group element e lA as the product e lAH e lAo °, where Ah is an element of the 
Heisenberg algebra, and Aqo = J2 n >i ^ (here we suppress the indices ji, j n ). The 

gauge transformation of Ck with the gauge group element e tAo ° leaves the component C^ 
untouched so it does not violate the required condition C^ = a k - Thus, modulo gauge 
transformations, we may write C& = — ie lQ,H dke~ lQ,H , where fin is an element from the 
Heisenberg algebra, or, equivalently, 

C k = OL k + A k = a k (q) + akjiq)^ , a kj u Jl a in = u kn , fya w = d k aji . (12) 

The two latter relations follow from Eq.([ll]). If Qi = Q^(q) are Darboux coordinates for 
the symplectic structure in question, then in local coordinates = djQ n 00 nk - 

The gauge arbitrariness of the solution (^) is not yet exhausted. The gauge transfor- 
mations from the Heisenberg group (A = A#) still apply. Representing the gauge group 

element as the product of operators e l j (no summation over j) and the phase factor 
e jA(0> (associated with the unit operator in the Heisenberg algebra), we calculate the gauge 
transformation of the solution (|12"D explicitly. Assuming some set of the Darboux coor- 
dinates to represent the functions aj k , we get that Qi —+ Q\ = Qi+ uj under the 
gauge transformation. The condition that C k = «fe should not be violated under any 
gauge transformation would yield that dkX^ = ^(A^^Q) (we omit the details) which 
specifies \(°\ The integrability condition for this equation leads to 

diQ\ uj kidjQx = uJij . (13) 

Thus, Q 3 X is just another set of Darboux variables. This, of course, might be expected at 
the very beginning since the residual gauge transformations do not change the structure 
of the solution flT2"|). 



The gauge freedom in choosing the Darboux variables appears desirable. Whatever set 
of Darboux variables we use in practical calculations of the physical amplitudes, the result 
will not depend on it. Changing from one set of Darboux coordinates to another implies 
a corresponding modification of the gauge invariance condition. The physical states and 
operators are modified so that the physical amplitudes remain unchanged. This freedom 
also allows one to make the quantum theory globally well-defined 0. 

4. The invariance of the theory under canonical transformations of the Darboux 
variables is crucial for the path integral formalism because one has to assume a specific 
form of the constraints in order to calculate the path integral. So we pick out some set 
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of Darboux coordinates QK The Hamiltonian of the effective gauge theory satisfying the 
conversion equation fllCf ) has the form 

H g = H0) , # = sym<f (0 + Q{q)) , (14) 

where the functions q % = q l {Q) are inverse to Q l (q), i.e., q l {Q(q)) = q % ■ We also assume 
some appropriate ordering of the operators d l to make the Hamiltonian hermitian. Making 
use of the explicit form of the constraints, it is easy to prove that [frj, d k ] = 0. The required 
boundary condition follows from the relation fi^ 1 = q\ 

The commutation relations of all canonical variables in the effective Abelian gauge 
theory are the standard ones, so we can construct the Wiener measure regularized path 
integral for the transition amplitude on the extended phase space along the lines discussed 
in Section 1. To simplify notation, we denote by b a the set q l ,Pi and 8 l , and will write 
{b a , b 13 } =uj a/3 . The amplitude K"(b, b') satisfies Eq.@, where q l are now replaced by b a 
and h(b) is the symbol of the operator ( |T4] ) in the canonical coherent state representation. 
The integral kernel K^(b,b') is not physical because it is not gauge invariant. To obtain 
the physical evolution operator, we make use of the projection method [§, ||. Consider 
a unitary operator U(£) = e^ J<Tj with £ J being parameters. For any physical state we 
have U(£)\i[)) = We define a projection operator V = J cIk(£)U(£) where the measure 
c?/c(£) is normalized to unity, / d/c(£) = 1. Then P'—P and V 2 = V. In general, the 
zero eigenvalue of the constraint operators may lie in the continuum spectrum of some or 
all the constraints. In this case we take a sequence of the rescaled projection operators 
csVs,cs > 0, as 5 — > 0, where V$ = f dKg(£)U(£), and the measure is chosen so that Vs 
projects the Hilbert space into a subspace where <5f < ^ 2 for < 8 « 1 0. 

So, our task is to develop a coordinate-free path integral representation of the gauge 
invariant transition amplitude (q\VK t V\q') = (q\KtP\q'). By construction, [H, = 0. 
To this end we consider the Schrodinger-like equation that depend on a set of general 
functions (, k (t) 

idtKtt® = [-j(b-B) 2 + h(B) + e(t)a k (B)^K^] (15) 

= {H v + e(t)* k (B)}Kttt} ■ 

As before, we have extended the original Hilbert space to an irreducible representation 
space of the canonical operators [B a , Pp] = iSp, [B a , B 13 } = [P a , Pp] = 0, while the original 

canonical operators now have a reducible representation b a = — u al3 Pp + ~B a . The 
functions h and a are the (lower) symbols of the Hamiltonian H and constraints &t in 
the original Hilbert space Tioo. With this choice we can prove the following important 
property of the solution to Eq. (|15"D with the intitial condition i^^ [^] = \ v 

(b\k v t \i\\b') = (blTe-^o^+^i^MM^ _^ (jb\Te- i fi dTl6+ * k <- T) * k] \U) (16) 

as v — > oo, for all t > 0. Here B k \b) = b k \b). Averaging Eq.(|i~6D over the functions £ fc we 
obtain the sought-for representation of the physical transition amplitude 

(b\K t V\b>) = Jim / dp(0 (b\kmW) , (17) 
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where dp(£) stands for a normalized measure for functions £ (r) on the interval r G 
[0, £]. In particular we may take the white-noise measure dp(£) = dp u (^) = Vt;e~^ J ^ dT , 
v — > oo, or a non-pinned Wiener measure As a point of fact, the only condition the 



measure should satisfy is that it must provide at least one action of the projection operator 
V within the time interval [0, £]. By going over to the path integral representation in the 
l.-h.s. of Eq. ([T7|), we derive our final result 

b(t)=b , 
6(0) =6' 

The path integral (0) with the measure dp(£) included has a genuine (finite and 
countably additive) measure which is covariant under general coordinate transformations. 
It provides the sought for path integral formalism for general Hamiltonian systems, which 
was our main goal. 

To relate the path integral ([TJJ) to the original phase space manifold, one can make a 
change of variables (which, in this case, is a well-defined procedure!) (q,p,0) — > (q, a, ■&) 
where d k = q k {9 + Q{q)) or, equivalently, 9 k = dq n a n i uj lk = J q dQ k for any integration 
contour. In this representation, the canonical pairs (q, a) are nonphysical degrees of 
freedom due to the translational gauge symmetry q — > q + 8£. The variables i? are gauge 
invariant and represent the physical degrees of freedom, and their Poisson bracket gives 
the original symplectic structure = ct/ 3 '(i?) and h = h($) by construction. If the 

integral over q, a and £ is done, the result will be given by the exponential of the classical 
action averaged with the induced non-canonical, non-flat Wiener measure on the original 
symplectic manifold. The induced measure will specify the unit operator kernel {dl'd') 
and generalized coherent-state representation of the non-canonical symplectic structure. 
This technical problem will be studied elsewhere. 

5. In conclusion, we would like to emphasize the following. A representation of the 
quantum mechanical evolution as a sum over paths should not depend on coordinates 
because the geometrical notion of a path on the symplectic manifold is coordinate inde- 
pendent (a path is specified by its tangent vector). Thus, in order to determine a formula- 
tion of quantum mechanics via the path integral on the symplectic manifold, or, in other 
words, to quantize a general Hamiltonian system on the manifold by the path integral 
method, it is fundamental that the corresponding functional measure is coordinate- free. 
The latter, in our opinion, requires an (induced) metric on the symplectic manifold. We 
have given a constructive way of how this regularization can be obtained. The procedure 
developed here provides a solution to the long-standing problem of constructing a rigor- 
ous coordinate-free path integral formalism on symplectic manifolds. The above idea can 
also be successfully applied to the path integral quantization of gauge and constrained 
systems, in general (in a flat phase space) |L0], 11]. The operator technique established in 



Section 4 can be used to greatly simplify the Wiener measure regularization of the path 
integral in gauge theories and second-class constrained systems proposed in our earlier 
works [ID, TT | . It is also allows one to generalize the result of [TDJ to most general (open) 



first-class constraint algebras. 
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Thus, the coordinate-free path integral formalism based on the Wiener measure regu- 
larization can be regarded as a universal, self-consistent, path-integral quantization tech- 
nique. 
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